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ABSTRACT 

We explore the non-BPS analog of 'AdS bubbles', which are regular spherically symmetric 
1/2 BPS geometries in type JIB supergravity. They have regular horizons and can be thought of 
as bubbling generalizations of non-extremal AdS black hole solutions in five-dimensional gauged 
supergravity. Due to the appearance of the Heun equation even at the linearized level, various 
approximation and numerical methods are needed in order to extract information about this system. 
We study how the vacuum expectation value and mass of a particular dimension two chiral primary 
operator depend on the temperature and chemical potential of the thermal Yang-Mills theory. In 
addition, the mass of the bubbling AdS black holes is computed. As is shown numerically, there are 
also non-BPS solitonic bubbles which are completely regular and arise from continuous deformations 
of BPS AdS bubbles. 
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1 Introduction 

Black holes in five dimensions have been extensively studied in the framework of both ungauged 
and gauged supergravity theories. Furthermore, many of these explicit studies have been performed 
in the context of the STU model, which corresponds to N = 2 supergravity coupled to two vector 
multiplets. Because of the presence of the graviphoton and two vector fields, the STU model 
generally admits three-charge black holes with up to two rotation parameters. In fact, in the 
ungauged context, many solutions have been constructed, and their explicit forms have often been 
proven to be useful, especially in the context of black hole thermodynamics and stringy microstate 
counting. 

At the same time, AdS black holes (and their variants) in gauged supergravity theories have 
found widespread application in the study of the AdS/CFT correspondence. For instance, an 
i?-charged black hole in global AdSs geometry corresponds to equilibrium non-zero temperature 

= 4 SU(A^) supersymmetric Yang-Mills theory on 5^ x R with finite chemical potential. Five- 
dimensional BPS 'black holes' were first constructed in [1], although there it was also realized that 
they are in fact naked singularities. Unlike the ungauged case, angular momentum must be turned 
on in order to obtain true BPS black holes in gauged supergravity. This was done in [3,4] for one 
rotation parameter, and subsequently generalized to two independent rotations. As in the ungauged 
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case, these black holes admit non-extremal generalizations [2], which can be further generalized to 
include rotations [5-9]. 

BPS objects play an important role in AdS/CFT, regardless of their precise nature, since they 
necessarily survive in both the strong and weak coupling regimes of the duality. Along these lines, 
the BPS naked singularities were shown in [10] to correspond to a distribution of giant gravitons in 
AdSs X S^, where they were also denoted 'superstars.' In terms of the Lin, Lunin and Maldaccna 
(LLM) boundary conditions for 1/2 BPS configurations [11], these superstars correspond to a disk 
in the phase space of free fermions with a uniform shade of gray, which may be interpreted as a 
coarse graining of an underlying spacetime foam picture of gravity [13]. 

The investigation of smooth 1/2 BPS geometries in [11] also led to the construction of a new 
set of 1/2 BPS 'AdS bubbles' which eliminate the naked singularity of the singular black holes- 
not through angular momentum but rather by turning on additional scalar fields ipi. These scalars 
are present in the full Af = 8 gauged supergravity but lie outside of the J\f = 2 truncation [11, 14]. 
These AdS bubbles are in fact completely regular geometries without horizons and, as such, admit 
a description in the LLM language as a deformation of the AdS disk into an ellipse in the phase 
space of free fermions. 

In order to generalize the above five-dimensional single-charge AdS bubbles to include three 
charges, the authors of [14] considered a consistent truncation of five-dimensional J\f = 8 gauged 
supergravity that retains the three U(l) gauge fields in the maximal torus of the SO (6) gauge group 
along with five scalar fields. Although this truncated system is not itself the bosonic sector of any 
supersymmetric theory, it nevertheless allows the construction and lifting of BPS solutions to yield 
supersymmetric configurations of the full M = 8 theory and hence the original IIB supergravity as 
well. 

Although LLM configurations are by their very nature pure states, and not thermal ones, it is 
natural to investigate the effect of turning on non-zero temperature starting from a particular LLM 
background. The main goal of this paper is to initiate such an exploration by constructing AdS 
bubbles away from the BPS limit. In particular, we will look for non-extremal AdS black holes 
with regular horizon and carrying (pi scalar deformations of the same sort encountered in the BPS 
bubbles of [11,14]. We will refer to these solutions as 'bubbling AdS black holes,' which are a subset 
of thermal AdS bubbles. Though not the focus of our paper, we will demonstrate numerically that 
there are also non-BPS solitonic AdS bubbles which are completely regular and horizon-free. 

We note that a coarse-graining of LLM geometries was considered in [12], and expanded upon 
in [13, 15], where the free fermion configuration was given an equilibrium non-zero temperature. 
Using the LLM correspondence, this configuration can then be mapped into a 'hyperstar' super- 
gravity background by appropriate transformation of the Fermi-Dirac distribution into the LLM 
z{xi,X2,y) function. It is not clear, however, whether this background actually describes a non- 
zero temperature field theory; it is by no means obvious that simply giving the matrix eigenvalues 
(the 'free fermions') a Fermi-Dirac profile would correspond to turning on non-zero temperature 
in A/" = 4 super Yang-Mills. If this were the case, then one would expect the corresponding su- 
pergravity background to have a horizon and to break supersymmetry. In [12], it was pointed out 
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that the coarse-graining has been taken over only the half-BPS sector of the full Hilbert space of 
type IIB supergravity, and that the result of using the complete Hilbert space should be a non-BPS 
background that has a horizon. Because of these considerations, we shall not pursue this direction 
of coarse-graining, but will look directly for bubbling AdS black holes in the supergravity itself. 

For the case of the AdS black hole, the main effect of moving away from the extremal limit 
is to introduce a 'renormalization' of the charge. One might hope that the AdS bubble might 
be generalized away from the BPS limit in a similar way, in which much of the original structure 
survives in a 'renormalized' form. Unfortunately, upon closer inspection of the equations of motion, 
it does not appear that the AdS bubble can be made non-extremal in such a simple manner. Thus, 
we must content ourselves with either approximate or numerical methods. 

We will proceed by performing a linearized analysis of the second-order equations of motion. 
Although the non-extremal AdS bubble solution is not known, explicit solutions are known in the 
two separate limits of either turning off non-extremality or turning off the (pi deformation. Since 
we will focus on solutions for which there is a horizon, we choose to consider the initial background 
to be that of the non-cxtrcmal i?-charged black hole. Wc will turn on the additional scalars ipi at 
linear order, which will then backrcact onto the other fields at higher order. 

The linearized (pi equations are second-order equations with four regular singular points, and 
hence may be mapped to the Heun equation. Unfortunately, in contrast with the hyper geometric 
equation which has three regular singular points, much less is known about the solutions to the 
Heun equation. In particular, the general two-point connection problem, i.e. relating local solutions 
at two regular singular points, remains unsolved. 

This is of course not the first time where the Heun equation has arisen in the study of the 
wave equation in AdS black hole backgrounds. In particular, calculations involving quasi-normal 
modes for black holes have generally led to expressions related to the Heun equation. In such 
cases, various methods have been applied in order to obtain approximate solutions of the resulting 
system. A common method, also used in black hole absorption calculations, involves matching 
approximate solutions in two overlapping regions: an inner region containing the horizon and an 
outer region that includes asymptotic infinity. Provided that there is a large overlap, then essentially 
complete information may be reliably extracted from this approach. Higher-order corrections can 
be included via a perturbative approach. We will implement this approximation technique for the 
regime T S> /Xj, where T and /Xj are the temperature and chemical potentials of the field theory, 
respectively. Furthermore, there is a high temperature limit T S> 1, along with T S> /Xj, for which 
there is an approximate solution which covers the entire region from the horizon to asymptotic 
infinity without the need for matching. 

The asymptotic behavior of (pi is related to perturbations away from the UV superconformal 
fixed point of the dual field theory. In particular, as in the BPS case, the normalizable mode of 
•pi corresponds to giving a vacuum expectation value (vev) to the dimension two chiral primary 
operator Zf, where Zi = -^{4>\ + '^4'2) ^| and ^2 three pairs of real scalars of = 4 super 
Yang- Mills theory, for i = 1, 2, 3. On the other hand, the non-normalizable mode of ipi corresponds 
to a term of the form Tr Zf in the Lagrangian of the conformal field theory. This massive term is 



4 



only present for nonzero temperature, which could be indicative of a phase transition in the field 
theory at zero temperature. We will discuss how these field theory deformations depend on the 
physical parameters of the thermal field theory, namely the temperature, chemical potential and 
i?-charge, in both the grand canonical ensemble and the canonical ensemble. 

The paper is organized as follows. In section 2, we review the previously-known AdS black hole 
and BPS bubble solutions. In section 3, we perform a linearized analysis of the bubbling AdS black 
hole. In particular, we consider the linearized ipi equations in the background of the AdS black 
hole. Focusing on the single-charge case, matching and perturbation techniques are used to find 
approximate solutions to the linearized (p\ equation. Properties of the (pi solution are discussed. 
Specifically, the asymptotic behavior is matched with perturbations of the dual field theory as 
functions of temperature and chemical potential. In section 4, we consider the backreaction of ipi 
onto the metric and other matter fields, and in section 5 we then discuss the mass of these non-BPS 
bubbles. In section 6, we present some numerical support for the existence of bubbling AdS black 
holes for arbitrary values of ipi . Lastly, we conclude in section 7. 

2 Review of AdS black holes and BPS bubbles 
2.1 AdS black hole 

The bosonic Lagrangian for the STU model takes the form 

2 3 

C = R*1-\Y^ *d(t)a A d^a - ^^Xr^*F' AF' - V*l + A A A^, (2.1) 

a=l i=l 

where the Xj's are constrained scalars satisfying XiX2X-i = 1, and which may be taken to be 

Xi = X2 = e-73^i+75^2^ = e^^\ (2.2) 

The scalar potential is given by 

3 

V = -Ag^ Yl = - V(^2X3 + X3X1 + X1X2), (2.3) 

i=l 

and the ungauged system is recovered by setting g = 0. 
The three-charge AdSs black hole solution is given by 

^ji) = -cothAif-^dt, Xi = n^/'^Hr\ <pi = 0, 

f = l-^+^VW, n = HiH2Hs, (2.4) 

where the harmonic functions are given by 

Hi = l + ^, qi = msinh^^i. (2.5) 
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We shall focus primarily on the case of a single charge qi, for which the roots of / are given by 

rl = -2^(1 + g^i) ± 2^ V(l + g^Qi? + ^a^rn . (2.6) 

The event horizon is located at = r+. Notice that, in the one-charge case with gi > 0, any 
positive value of m guarantees a regular horizon. On the other hand, for the two and three-charge 
regular horizon may always be obtained for sufficiently large m. 
This supergravity background is dual to equilibrium non-zero temperature Af = 4 SU(A'") su- 
persymmetric Yang-Mills theory on 5^ x M with chemical potentials for the U(l) i?-charges. The 
temperature of the field theory is equated with the Hawking temperature of the black hole, which 
is [20] 

^_ 2rl + ri{l + Z^q^)-n^1^ ^ ^2.7) 



where the horizon radius Vh is the largest root of /. Wc have set g = 1 for simplicity. Likewise, the 
i?-chargc chemical potentials //j of the field theory are equated with the electric potentials at the 
horizon, which are 

where 

= Qiirl + qi) + \ \l{rl + q^)] . (2.9) 

Also, the physical charges Qi of the AdS black hole correspond to i2-chargcs in the dual field theory. 
These relations will be useful for expressing various results in terms of the physical quantities of 
the field theory. 

2.2 BPS AdS bubble 

The authors of [14] considered a consistent truncation of five-dimensional Af = 8 gauged super- 
gravity retaining the three U(l) gauge fields in the maximal torus of the S0(6) gauge group along 
with five scalar fields. In A/" = 2 language, this corresponds to taking the bosonic sector of the STU 
model and coupling it to three additional scalars ipi, which are not described by special geometry. 
The Lagrangian is given by [14] 



a=l 1=1 
3 3 

i ^ *d(pi A difi - 2g'^ ^ sinh^ ipi*A^ A A\ (2.10) 



2 

i=l 1=1 

where the scalar potential has the modified form 

3 3 



V = 2g'^ Xf sinh^ (^j - 2 ^ XiXj cosh ipi cosh ipj^ . (2.11) 



1=1 i<j 
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The original STU model scalars satisfy the same constraints as above, and in particular may be 
given by 

In order to ensure a supersymmetric solution, the regular three-charge AdS bubble configuration 
then takes the form [14] 

ds^ = -n-^/^fdf +n^/^{f-^dr^ + r^dnl), 
vlji) = -H-^dt, Xi=n^/^Hr\ cosh^i = {xHiY, 

f = l + g^xH, n = H,H2Hs, (2.12) 

where x = r^, and a prime denotes a derivative with respect to x. The above solution is fully 
determined up to the functions Hi, which must satisfy the conditions 

fixHi)" = -g^xHi^ - l]nHr\ (2.13) 

to ensure that the equations of motion are satisfied. 

For the one-charge case, corresponding to a 1/2 BPS configuration, we may take H2 = = 1. 
Then the equation of motion (|2.13|) reduces to 

[{l + g^xH,fr = 2g\ (2.14) 

which admits a general solution of the form 



y g'^x g^x"^ g'^x 

Here ^1 is the i?-charge, and c is a constant related to the ipi scalar deformation. Note that, at 
large distances. Hi admits the expansion 

~ 1 + — + ••• , (2.16) 

X 

while regularity of the AdS bubble at short distances demands c = 1. In addition, the BPS naked 
singularity of [1,2] is recovered by taking c = 1 + g^qi, in which case Hi reduces to the standard 
'harmonic function' form Hi = 1 + qi/x. 

In general, as we will demonstrate below, the scalar ipi carries Eq = 2, and hence has an 
asymptotic expansion of the form 

Cl + C2 log X 

<^i ~ \ , (2.17) 

X 

However, based on the explicit solution ()2.15p . we see that the log term vanishes, C2 = 0, while 



Ci = ^{l+qif-7? = yJqi{qi + 2), (2.18) 

where the second expression is for the regular AdS bubble solution. We have set g = 1- Applying 
the AdS/CFT dictionary (see e.g. [19]), this implies that the dimension two chiral primary operator 
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Tr Z\ gets a vev ci, where Z\ = "^('/'i +i'/*2) and (f)i and (/)2 are two of the six real scalars of = 4 
super Yang-Mihs theory. 

Although we are not aware of any closed form expressions, it can be seen that regular two and 
three-charge solutions to (|2.13p also exist. These solutions correspond to 1/4 and 1/8 BPS AdS 
bubbles, and can be realized in the framework of LLM configurations with fewer supersymmetries 
[16-18]. 



3 Linearized analysis of the bubbling AdS black hole 



Given the relatively simple non-extremal generalization (j2.4p of the BPS AdS black hole solution, we 
have been led to look for a corresponding non-extremal version of the AdS bubble (|2.12|) where the 

scalars are present. In the absence of supersymmetry, we no longer have the benefit of working 
with first-order Killing spinor equations. However, the form of the non-extremal black hole (j2.4p 
is curiously close to that of the BPS limit; the primary difference is that the charge parameters 
qi in the harmonic functions Hi are 'renormalized' as follows: qi msinh^/3i. This suggests that 
perhaps a similarly straightforward generalization may be obtained for the AdS bubble, where much 
of the BPS structure might survive, except perhaps in 'renormalized' form. 

Unfortunately, closer inspection of the equations of motion arising from the Lagrangian (j2.10p 
does not suggest any simple manner in which the AdS bubble solution may be made non-extremal. 
In particular, retaining the BPS-like relation cosh(/7j = {xHi)' leads to either the possibility that 
'^i = 0, in which case Hi takes the form (j2.5p . or to a constrained system of equations which only 
appear to admit a natural solution of the BPS bubble form given by (|2.13p . As a result, we must 
use either approximate or numerical methods when moving away from extremality. 

3.1 Linearized (pi equations 

In this section, we explore the basic features of the non-extremal AdS bubble by performing a 
linearized analysis of the second-order equations of motion obtained from (j2.1Up . We proceed by 
noting that, although it is not clear how to write down the complete solution for a non-extremal 
AdS bubble, explicit forms for the solutions are known in the two separate limits of either turning 
off non-extremality or turning off the ipi deformation. Since we are mainly interested in solutions 
with a horizon, we choose to start from the non-extremal i?-charged black hole background of (|2.4p 
and (|2.5p and then turn on the additional scalars ipi at linear order. Turning on these scalars will 
then backreact onto the other fields. However this backreaction occurs at the next order, and may 
be ignored in the initial analysis. (We will return to the backreaction in another section.) 

At linearized order we are only concerned with the linearized equation of motion for which 
takes the form 

3 

□ - Ag''Xi{2Xi -Y^Xj) - 452(4)2] = 0. (3.1) 



Here the scalars Xi take on the background values given above in ()2.4p and (j2.5p . There is a 
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slight subtlety for the background gauge fields, however, related to the 'gauge fixed' form of the 
action (j2.10p . In particular, while the standard black hole solution ()2.4p is invariant under gauge 
transformations of A^, the above equation of motion is not. Taking this into account, we allow a 
constant term in the expression for the background electric potential 



A' = {b' - coth piHr^)dt = Ai dt. 



(3.2) 



We shall see below that this constant must be chosen to make the potential vanish at the horizon. 
The motivation for this requirement can already be seen by noting that the invariant square of 
the electric potential, (vl^)^, which acts as a source in (|3.ip . blows up at the horizon unless A\ is 
arranged to vanish there. 

By substituting the background fields into (|3.ip . we obtain the scalar equation 



+ — - 2b'^yl + m/q^x + {b'fix + qi 



= 0,(3.3) 



where 



+ Qj) + 



mx. 



is a cubic polynomial, and where qi = msinh^ Note that, in the event qi 
is replaced by the considerably simpler expression 



(3.4) 

0, the above equation 



{x'm' + '-(x'-llq, 



0, 



(3.5) 



which may be obtained by directly taking A* = in (j3.ip . 

Since x^/ is cubic, it can be seen that (j3.3p is a second-order equation with three regular singular 
points at the roots of x'^f. Including the singular point at infinity, which is also regular, this 
equation in fact has four regular singular points, and hence may be mapped into the general Heun 
equation. Unfortunately, the general two-point connection problem is as yet unsolved. However, 
general features of the linearized (pi deformation may be extracted from the second-order equation 

(1231). 

Recalling that the goal of the linearized analysis is to turn on a ipi deformation starting from the 
i?-charged black hole background, we demand that the solution to (j3.3p be regular and bounded 
in the entire region from the horizon to the spatial boundary at infinity. Before examining the 
solution, we find it convenient to trade the non-extremality parameter m with the horizon location 
Xh, defined to be the largest positive root of /(x). Note that, in terms of x/j, we have 



m 



Xh 

1=1 



(3.6) 
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along with the factorization 

x^f = g\x-Xh){x^ + Qx-R), (3.7) 

where 

Q = a;^ + ^ + + 92 + 93, -K= • (3.8) 

Examination of the indicial equation around the horizon (which is a regular singular point) 
yields the characteristic exponents ±Ci where 

Noting that the expression above is non-positive (since Xh + % must be positive to avoid naked 
singularities) , we immediately see that the characteristic exponents are purely imaginary, except for 
the case when they vanish. Since imaginary exponents give rise to undesirable oscillatory solutions 
of the form Lpi ~ sin(|^j| log(x — j;^)), we conclude that the constant 5* in the electric potential must 
be adjusted to satisfy 

V qiXh + qi Hi{xh) 

This ensures that the potentials given by (|3.2p indeed vanish at the horizon, A^{xh) = 0, thus 
confirming what we had alluded to above. 

Demanding that the electric potentials vanish at the horizon, the ^pi equation can now be put 
into the form 



{{x - Xh){x^ + Qx - R)^[)' + 
'xl + Qxh - R 



^{x + qj) - {x^ + Qx - R) ^Yl{x + qj] 



+ (X - Xh] 



^i = 0. (3.11) 



Xh + Qi 

The solutions to this equation may be characterized by the Riemann P-symbol 



Xh X\ X2 OO 



P < 



ai 02 1 ; 2^ y ! (3.12) 
—ai —a2 1 



where xi and X2 are the two roots of the quadratic equation + Qx — R = 0. Their associated 
exponents, ai and 02, may easily be obtained, although we have no particular need for their explicit 
forms. We do note, however, that the characteristic exponents sum to 2, which is in agreement 
with the general theory of the Heun equation. 

The single-charge case 

Since we are primarily interested in the non-extremal generalization of the single-charge 1/2 
BPS bubble, we now focus on the case where 92 = 93 = 0. In this case, we assume the corresponding 
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scalars ip2 and 933 ought not to be excited. As a result, we are left with a single equation for 
which now takes the form 

Xh + Q qi 



{x{x - Xh){x + Q)(pi) +x 



1 + 



fi = 0. (3.13) 



x + Q Xh + qi_ 

This equation still contains four regular singular points, and as such can be brought into Heun 
form. In particular, we may map the singular points {0,Xh, —Q} to {0, l,a} by introducing 

X = zxh, Q = —axfi, (3-14) 

after which the scalar equation becomes 



{z{z - l){z - a)ip[y + z 



ipi = 0. (3.15) 



z — a 

Here the prime denotes a derivative with respect to z. We have also defined the dimensionless 
parameter 

C= . (3.16) 

Note that ^ — > 1 in the limit of large charge, while ^ ^ in the limit of vanishing charge. 

The scalar equation ()3.15p can be brought into the canonical form of the Heun equation through 
the substitution ipi = {z — a)^^(fi. As a result, the solution may be written in terms of a local 
Heun function 

^1 = ^Hl{-, 1 - 1 + v^, 1, 1; ^). (3.17) 

z — a a a z — a 

Here we have imposed the boundary condition that (pi is regular at the horizon. In particular, the 

expansion of the local Heun function gives 

1 I n I r\'-^ I 4-a + (6-a)C + CM^-i)^ , ..... 

'/'i = l + (l + C)^ + 4 7^31)2 +•••• (3-18) 



where C, is given in ()3.16p . and where 



1 + ^1±^~). (3.19) 

Xh 



In principle, we would also like to obtain the expansion of ipi in the asymptotic regime 2; — > c« in 
order to extract the boundary behavior (j3.20p . However, the two-point connection problem for the 
Heun equation is in general a difficult task, and there is as yet no straightforward way to connect 
the behavior of (j3.17p near the horizon with (j3.20p near the boundary. In order to match the 
horizon and boundary behaviors of we need to utilize further approximation techniques. Before 
turning to this, we shall first consider the general behavior of 934. 

General behavior 

The asymptotic behavior of the solution for (pi is governed by the point at infinity. In particular, 
given the repeated characteristic exponent of 1, it may be shown that the solution near the boundary 
has the form 

c\+c\\ogx 

(pi{x) ~ \ , X ^ 00, (3.20) 
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which is consistent with (pi being an Eq = 2 scalar (or equivalently being associated with conformal 
dimension A = 2 in the dual gauge theory). Although the two independent solutions encoded in 
(j3.20|) admit different interpretations in the dual gauge theory, they are both allowed to be present 
in the supergravity solution. 

In particular, as in BPS case, c\ can be interpreted as the vev of the dimension 2 chiral primary 
operator Tr Zf, where Zi = -^{4>i + ^4>2) 4'\ ^^'^ 4'2 three pairs of real scalars of = 4 
super Yang-Mills theory, for i = 1,2,3. On the other hand, the c| term corresponds to adding 
the relevant deformation c| Tr Zf to the Lagrangian of the conformal field theory. The term is 
only present when one moves away from extremality, which indicates that there may be a phase 
transition in the field theory at zero temperature. 

Turning to the horizon, we see that, with the choice of 6* in (|3.1Up . the repeated characteristic 
exponent of at x = x/i indicates that the solution near the horizon has the form 

V9j(x) ~ + 41og(x - x/j) H , x^Xh- (3.21) 

Clearly, we must set d\ = to avoid a logarithmic divergence of Lfi at the horizon. 

In principle, this boundary condition now fixes the complete solution, in the sense that the 
coefficients c\ and C2 in the asymptotic expansion may be determined directly from d\ and the 
physical parameters of the background solution. However, in the absence of any general connection 
formulae for the Heun equation, there is no straightforward way to make this relation explicit. (We 
note that, even if the generic connection matrix were known, it might not be applicable to this 
solution as it has repeated characteristic exponents.) In order to match the horizon and boundary 
behaviors of Lpi for the single-charge case, we instead turn to approximate solutions of equation 
(I3T5]) . 

3.2 Matching approximate solutions 

We first consider an approximation technique which involves matching approximate solutions in 
two overlapping regions. In particular, one region contains the horizon while the other includes 
asymptotic infinity. Provided these two regions overlap, the solutions may then be matched up 
in the overlap region. A drawback of this approach is that it is not always possible to ensure a 
large overlap region, depending on the physical parameters of the system. However, if such a large 
overlap exists, then essentially complete information may be reliably extracted using this matching. 

Noting that a introduces a new scale into the problem, we may consider solving (j3.15p in the 
two regions: i) the asymptotic region where z ^ 1, and ii) the horizon region where z <C \a\. 
Provided \a\ ^> 1, these two regions will have a large overlap (1 <^ z <^ \a\) where reliable matching 
may be performed. From (|3.19p . we see that overlap is ensured for either qi ^ Xh or x/i <C 1/5^. 

3.2.1 The asymptotic region 

To highlight the asymptotic region, 2; ^ 1, we may rewrite the scalar equation (j3.15p as 
\w w — 1 J w{w — 1)V w — I J a lw{w — 1/a) \ w — I ) \ 
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where we have introduced the rescaled coordinate w = z/a. In this case, the additional factor of 
1/a on the right-hand side of (j3.22|) allows us formally to develop a solution for ipi as a perturbative 
expansion in 1/a 

= ^(0) + i^(i) + J_^{2) . . . , (3.23) 
a 

where (p^'^^ solves the homogeneous equation corresponding to the left-hand side of (j3.22p . Since 
this can be put into hypergeometric form, the solution is essentially known. In practice, however, 
matching of the asymptotic and horizon expansions is facilitated by introducing yet another expan- 
sion, this time in C,. Examination of (j3.16p indicates that there are two relevant limits to consider, 
namely the large and small charge limits. 

The large charge limit corresponds to qi » Xh, or equivalently ^ w 1. In this case, we let 
C = 1 + C and rearrange (I3.22[) to read 



\w w — 1 J 



2 I \ , 'w-2 C 1 



fl + Tf^ 



■w{w — 1/a) I w — 1 



w{w — 1)^ w{w — 1)2 a 

(3.24) 

Solutions to this equation can now be developed as a double expansion in and 1/a. Although 
this may seem to be only a slight rearrangement of (j3.22p . the main simplification here is that 
the homogeneous equation can now be solved in terms of elementary functions, ^p^^^ = ciui -|- C2M2 
where 

ui = U2 = (log{-w) + - ) . (3.25) 

w — 1 w — 1 \ w J 

At each successive order in the perturbation, the lower order solutions feed in as sources on the right- 
hand side of p.24p . However, since the homogeneous solutions are elementary, the inhomogeneous 
system has a straightforward solution which can be developed, e.g through variation of parameters. 

Up to first order in both ( and 1/a, we find that the two linearly independent solutions for the 
outside function can be expressed as 



V'out 



2 

V>out 



1 



W — 1 
1 

■w — 1 



log{-w) + - - |flog(-u;)(log(l - 1) - 1) + (2 - i) log(l - - 

+ log2 (1 - -) + 2Li2 (-^) log{-w) + + Oio) + 
\ wy \1 — wy / aw \ 2w / 



(3.26) 



where we recall that w = z/a. We have organized these two solutions according to the asymptotic 
behavior 

^out~^ + ---, </'out~^logz + --- , (3.27) 
corresponding to the normalizable and non-normalizable Eq = 2 modes of (pi, as in ()3.20p . 



13 



3.2.2 The horizon region 

Turning next to the horizon region, z <C \a\, we now choose to write the scalar equation ()3.15p as 



z z — 1 



1 



1 



^1 + 



1 



z- 1 



I — z/a 



(3.28) 



a \1 — z/a^ 

This again allows us to develop an expansion in 1/a, where the independent solutions to the 
homogeneous equation on the left-hand side are simply 

z - 1 



Ul = 1, 



U2 = log 



(3.29) 



These two solutions correspond directly to the near horizon behavior given in (j3.2ip : in particular, 
we see that only ui is well behaved at the horizon. Developing this solution to the first few orders 
in 1/a gives 

V^in = l + ^{l + C){z-'l + logz) + ^\{z-l){9{l-e) + 2{e + 7C + 4)z) 



24a2 



- (5C + 14C + 5 - 6(1 + C)'^) log z + 3(1 + C)' log' z + 6(1 + O'Lia [1 - - 



+ 



(3.30) 



Note that this is purely an expansion in 1/a, and in particular it is valid for arbitrary ^. Further- 
more, it can be seen that this expression for (^in agrees with the near horizon expansion given in 
(j3.18p in the overlapping region of validity z ^ 1 and a — > cxd. The advantage of (j3.30p over (j3.18p . 
however, is that here (pi^ remains valid even for z away from the horizon (provided z <^ \a\). This 
is precisely what is needed in order to match the horizon expression with the asymptotic forms of 
the solution given above in (|3.26p . 



3.2.3 Matching 



While the asymptotic and horizon solutions (|3.26p and ()3.30p were derived under the independent 
conditions of z ^ 1 and z ^ \a\, they are both valid in the overlap region 1 <C z <C |a|, so long 
as \a\ ^ 1. In order to match the solutions in this overlap region, it is convenient to rewrite the 
horizon solution (j3.30p in terms of w. 

After suitable rearrangement, (pi^ then takes the form of a series solution in w as well as an 
expansion in 1/a 



W 



1 + f (1 + C) + ^iC + 7C + 4) + 0{w-') 



1 

+- 
a 



w 



-(1 + C)(log(u'a) - 1) + —(-lie - 14C + 1 + 6(1 + C)' log(u'a)) + 0{w' 



24 



(3.31) 



Here it is important to realize that, although the expressions are no longer complete at each new 
order in 1/a, the resulting series in w are in principle well behaved for \w\ < 1. This is what allows 
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a consistent matching with the asymptotic solution (j3.26p . with (/9out similarly expanded as a series 
in w. 

The asymptotic solution to the scalar equation is in general a linear combination of the two 
solutions given in (|3.26p : 

/X 1,2 axh{ci+C2\og{-x/axh)) „„x 
V?out(a:) = ciV9out + C2V?out > (3-32) 

X 

where we have transformed back to the coordinate x = r^. By matching this with the normalized 
horizon solution (I3.3ip . we find 



i (l - log(-<.) - |(y + 1 + log(-»))) + 0(?) + 0(i 



where 



C = (3.34) 

and a is given by (j3.19p . Note that both ( and a are negative for physical values of the charge and 
horizon radius. 

In order to more readily apply the AdS/CFT dictionary for examining the asymptotic behavior, 
we express ()3.32p in terms of the original r coordinate: 

axh[ci - C2log{-axh)] logr 
^out{r) = ^ h 2axhC2 ■ (3.35) 

Note that we can always include factors of g to ensure that the logarithms have dimensionless 
arguments. However, we have set g = I for simplicity. The above term corresponds to giving 
a vev to the chiral primary operator Ti Z^. The (log r)/r^ term, on the other hand, corresponds 
to adding a massive deformation —^m^TrZ^ to the Lagrangian of the conformal field theory. A 
priori, there is an ambiguity in the normalization of (pi. Recall that, for the BPS AdS bubble, the 
thermal mass fh = while, from (|2T8]) . TrZ^ gets a vev vi = \J qxicjx + 2). Therefore, since we are 
interested in a thermalization of = 4 super- Yang Mills on top of the 1 /2 BPS sector specified by 
a given vev v\, we may fix the normalization by taking (TrZ^) = v\ for the bubbling AdS black 
hole. Then, taking the limit q\ ^ Xh, we find 

^ Axh - r^) y^qiiqi + 2) . (3.36) 

As Xh vanishes, the thermal mass goes to zero. This is certainly not unexpected, since the thermal 
mass vanishes in the BPS limit, though at the same time it would not have been surprising if in 
did not behave in a smooth manner in this limit. Interestingly enough, fh also vanishes for qi — > 1. 
As we will see, this could be due to a phase transition in the field theory that is associated with 
the Hawking-Page transition. In particular, for qi <1, the operator Tr could become tachyonic, 
thereby destabilizing the moduli space of the field theory. 
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We would now like to express m in terms of the physical parameters of the thermal field 
theory, namely the temperatm'e T, chemical potential /ii and ii-charge Qi. These quantities are 
all expressed in terms of the AdS black hole parameters qi and = in (|2.7|) , (j2.8|) and (|2.9|) . In 
the limit qi » Xh for a single charge, 

T = I , /ii = \/l + x/i , Qi = giVl + Xh . (3.37) 

We therefore have a choice of expressing Xh and gi in terms of two out of the three physical 
parameters of the field theory, which corresponds to different ensembles. For example, the -R-charge 
Qi is held fixed in the canonical ensemble. Therefore, the thermal mass is a function of temperature 
and chemical potential in the grand canonical ensemble, while it depends on temperature and R- 
charge in the canonical ensemble. 

We first work in the grand canonical ensemble. We can express Xh and qi in terms of T and /xi 
in the limit qi ^ x/^ as x/j « /if — 1 and qi ~ 2txT{ttT it Vtt^T^ — 1) — 1. Notice that black holes 
with two different values of qi correspond to the same temperature. However, since the entropy is 
given by S* ~ 2TTXh^/qi, the black hole with larger qi is entropically favorable and so we shall take 
the + sign in the expression for qi. Since Xh > 0, in this regime fif > 1. Also, vrT > 1 in order for 
qi to be real. 

One way in which qi ^ x^ is if T ^ fii. This implies that T ^ 1 and qi ~ 47r^T^. We can then 
express the thermal mass as ~ 2(/if — 1). Thus, for a given chemical potential, the thermal 
mass vanishes at zero temperature and approaches a constant in this high temperature limit, which 
may imply a sort of saturation taking place. Also, the thermal mass gets enhanced by increased 
chemical potential. Notice that this high temperature regime does not include the point at which 
the thermal mass vanishes at a finite temperature. 

Alternatively, qi ^ Xh can be satisfied without having to take T to be large by taking /Ui ~ 1. 
Then after expressing the thermal mass in (j3.36p in terms of T and fii, we find that fh = for 
ttT = 1, which corresponds to gi = 1. As T increases, the thermal mass asymptotically approaches 
the constant value discussed in the previous large temperature limit, that is 2(/if — 1). 

We will now look at the situation in the canonical ensemble, for which the -R-charge Qi is held 
fixed. One way to satisfy qi ^ x^ is to consider ^ Qi and T ^ 1. Then we have 

where > since x^ > 0. We see that the thermal mass decreases with temperature and increases 
with the -R-charge Qi. Note that in order for there to be a horizon at Xh > 0, fh? > 0. This requires 
that Qi > Att'^T'^. 

We will now briefly discuss the conditions for the regime qi S> Xh to be consistent with local 
thermodynamical stability constraints and dynamical considerations for the charged AdS black hole 
that were given in [20], in both the canonical ensemble and grand canonical ensemble. At the level 
of linearized ifi, we do not need to consider the backreaction on this background, which would then 
alter these thermodynamical relations. 
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In the grand canonical ensemble, in order for the AdS black hole to be dynamically preferred 
over pure AdS, vrT > vrTc ~ 1, where Tc is the temperature of the Hawking-Page transition. 
Furthermore, the local thermodynamical stability constraint is satisfied only for vrT ~ 1. Thus, 
we only have a small window for which the above calculations are consistent with stability. The 
temperature at which in vanishes lies within this window; in particular, this occurs at the Hawking- 
Page transition and presumably signifies the corresponding phase transition in the field theory. On 
the other hand, the large temperature regime where ifi saturates does not lie within this window. 
For the canonical ensemble, the local stability constraint is satisfied for vrT ~ ^/2, which also 
satisfies the condition for the AdS black hole to be dynamically preferred over pure AdS. 

We would now like to recall the approximations that have been made, namely that \C,\ <^ 1 and 
\a\ S> 1 correspond to the large charge limit qi ^ x^, which is equivalent to T ^ i^i. We could also 
match approximate solutions in the small charge limit qi <C x/j, which means that |C| *C 1. Then 
in order to be consistent with the condition that the asymptotic and horizon regions have a large 
overlap, we require that \a\ ^ 1, which further implies that x/i <C 1. We can use ()2.7p and (j2.8p 
to express these conditions as //i <C 1 <C T. Therefore, this is the regime of high temperature and 
small chemical potential. However, as we will see in the next section, we do not have to rely on the 
matching technique if T ^ 1 and T ^ ni, which encompasses the above high temperature regime. 

3.3 Perturbative approach for a second high temperature regime 

We will now consider the linearized ipi equation (|3.15p for the case of Xh ^ qi and Xh ^ I, where 
we have taken g = 1 for simplicity. As we will see, this is a second high temperature regime. Note 
that ipi can be expanded for small xj^^ as 

=uo(l + 1 ui + 0{xf )^ . (3.39) 

The function uq satisfies the following linear differential equation 

{{z'^ -l)zu'o)' + zuo = 0, (3.40) 

which has the general solution 

uo = — K{1 -z^) + ^ K{z^) . (3.41) 
vr vr 

where K{x) denotes the complete elliptic integral of the first kind. For the solution to be regular 
at the horizon z = 1, it is necessary to set C2 = 0. Furthermore, without loss of generality we can 
set the scaling factor ci = 1, so that near the horizon no behaves like 

no = 1 - i(z - 1) + ^(z - If + 0{{z - If) . (3.42) 

Asymptotically we find that 

2[log^-7-^(l/2)] 
no , (3.43) 

TTZ 

where 7 is Euler's constant and ip is the digamma function. Approximate numerical values are 
7 « .58 and '4>{l/2) ^ -1.96. 
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We can now consider the x^^ corrections in ipi. The function ni can be expressed in terms of 
a functional integral as 

r v{y) 

V = I+iz) + qil^iz) , ^±(^) = / :^{y{y-l)uo{y)4{y)±yuliy)yy. (3.44) 



Setting the lower bound of the v integral to unity ensures that the solution remains regular at the 
horizon. Moreover, the chosen lower bound of the ui integral guarantees that ui vanishes at z = 1, 
so that the boundary condition on ipi remains unchanged. 

To compute the asymptotic behavior in ui, evaluate /±(oo) in order to take the large z expansion 
of the Ul integral. The result gives 

7r(/+(oo) + gi/-(cx))) 

(fi ^uo , (3.45) 

2xhZ 

for large z. The approximate numerical values are I+(cxd) 1.03 and I-{oo) ~ —3.51. In terms of 
the coordinate r, the asymptotic behavior is roughly given by 

^1 - l^^ii^ + 4 , c ^ -^[logx, + 7 + ^^(1/2)] - 5[/+(oo) + (Zil-(oo)] . (3.46) 

Once again, the 1/r^ term corresponds to giving a vev to the chiral primary operator TrZ^. We 
normalize by equating this vev to the value vi = \J qi{qi + 2) for the case of the BPS AdS 
bubble, as we did in the previous section. This enables us to find the thermal mass of Tr to be 



^ V^M±1 , (3.47) 

log Xh 

in the limit Xh ^ qi and x/^ ^ 1. We would like to express this in terms of the temperature, 
chemical potential and i?-charge of the field theory. From (|2.7|) and (|2.8|) and (|2.9|) . we can express 
the temperature, chemical potential and i?-charge Qi in terms of the black hole parameters x^ and 
qi in the regime Xh ^ qi and x/i ^> 1 as 



Tss^^^, /UiPs^, Qi^y/qiXh- (3.48) 
vr 

In the grand canonical ensemble, this regime corresponds to T ^ ^ui and T ^ 1 and the thermal 
mass can be expressed as 

^,^V5!M±1. (3.49) 

2iogr ^ ^ 

Thus, the thermal mass increases with chemical potential and decreases with temperature. 
In the canonical ensemble, this regime corresponds to T'^ ^ Qi and T ^ 1 and 



2vr4r4iogr • ^"^-^"^ 

The thermal mass now increases with i?-charge Qi and still decreases with temperature. Note that, 
in this regime, the results for the thermal mass are identical for the case of three equal charges. 
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Recall that there was only a small window for which the regime discussed in the previous section 
satisfied the local thermodynamic stability constraints and entropic considerations. On the other 
hand, for the regime 3> qi and Xf^ ^ 1, all of these constraints are satisfied. This provides a 
large range of temperatures for which this system can be reliably discussed. Of course, backreaction 
and non-linear effects would also have to be taken into account, should we desire a more detailed 
treatment of the thermodynamics. 

4 Taking backreaction into account 

In the previous section, we have explored the linearized equations of motion for the deformation 
scalars ipi. Going beyond linear order, these scalars will backreact on the metric through the 
Einstein equation, as well as on the other matter fields through the couplings implicit in the 
Lagrangian (j2.10p . In particular, the field ipi acts as a source for the gauge fields through the 
sinh^ ipi * A A'' couplings, and the Xi (or equivalently (/>i and (f>2) scalars through the modified 
scalar potential (|2.1ip . 

While the generalized backreaction equations are straightforward to obtain, we restrict the 
analysis to the non-extremal generalization of the single-charge AdS bubble. In this case, the 
natural way to parameterize the metric backreaction is to start from the (one-charge) black hole 
solution (|2.4p and to write 

ds^ = -H-^/^fdt^ + H^/^{f-^dr^ + r^dnj), (4.1) 
where the metric functions are, to second order in the linearization parameter e, 
H = l + ^ + e2/,2(x) + ... , 

X 

f = l-!I^+g^xH + e^f2{x) + --- = l + g\i-- + g^x + eHf2+9^xh2) + --- . (4.2) 

X X 

Here we recall that x = and qi = m sinh^ Pi , as indicated in ()2.5p . Note that we have continued 
to write the function / in the natural combination of / = 1 — m/x + g^xH plus corrections. 
We now turn to the matter fields. For the electric potential, we follow (|3.2p and write 

Aj = {b^ -cothpiH-^) +e^a2(x) + ■■■ 
, n f Xh X \ 2 

= coth Pi ■ ■ ■ ^ e 02 H 

\xh + qi x + qi+xe^h2j 

= -gicoth/?!- -I- ( coth /?! r5-/i2 + 02 H , (4.3) 

{xh + qi)[x + qi) V [x + qiY J 

where coth f3i = y^l + m/qi = \J {g"^ /qi){xh + l/g'^){xh + qi)- Finally, for the backreaction of ipi 
on the Xi scalars, we note that it is consistent to set X2 = X^ = 1/yJXi, since ^pi sources X2 and 
X3 in an identical manner. We then take a multiplicative parameterization of Xi: 

X, = H-^l\l + e\2{x)) = (1 + 1)"^' (1 + (x2 - l^/^^)) + • • • • (4.4) 
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For Lp\ of order e, its backreaction on the metric fields /12 and /2, the electric potential 02 and 
scalar X2 is then governed by the set of inhomogeneous second-order equations 

{xh + + X/, + l/g^ + Qi)^ 
+ g^x^K + 2x2/i2]' + 5'[3^'/^'2 + 2<ZiX2] 

_ + X,, + 1/52 + 2gi + + l/g^)/{xh + qi)) 2 



X + Xh + 1/5-2 



(4.6) 



((x + (7i)^4)' = coth/?ix(x + qi) ( (^'2 _ ^ ^"^^ 1/ 2^ ^2 '^i) ' (4.7) 

V [xh + qi){x + Xh + l/g^ + qiY J 

[(x(x - Xh){x + Xh + 1/5^ + 9i)x2)' + 2^X2] + [2xf2/g^ + x^h2\' + [x^h'^ + 2giX2] = x^\. (4.8) 

In addition, the constraint equation coming from the Einstein equations gives rise to a rather 
cumbersome first-order equation 

■ X - Xh){x + Xh + I g + qi)X2 + ^g qi \ 7 \ yj ^ 

x + qi V [x + qiY J 

x{2,x + 2qi) ( qf n\ f ^ ^ a , 

\ /2 + 7 — \ NO - 3 /2 + 2gi coth/3ia2 

x + qi \{x + qiY J 

. 2 2 xxh{i + g'^{xh + qi)) \ ^1 ( 2 qixh{'^ + g'^{xh + qi)) 

45 X ; j ri2 — \ ^g X 



x + qi J V {x + qiY 

-g'^x{x - Xh){x + Xh + l/g^ + qi)v'i 
g'^x{x + Xh + l/g^ + 2qi + qi{xh + l/g^)/{xh + qi)) 2 



{x + Xh + l/g^ + qi] 



ft- (4.9) 

We observe that the metric backreaction equations (j4.5p and (j4.6p . along with the Xi equation 
(j4.8p . are coupled in a non-trivial manner. On the other hand, the electric potential equation (j4.7p 
is independent, and can be solved by quadratures: 

dx' , , „ 

a2(x) = coth/3i / — -TT / dx"x"(x" + gi) X 

J [x' + qiYJ 

V'hx")- , , , .N. ^f(^^OV (4.10) 



(x;, + gi)(x" + X;, + 1/(72 + ^1)2 ^ 

Note that the indefinite integrals allow for the addition of an arbitrary homogeneous solution 

a2 = ki + k2 — ^- — + ■■■. (4.11) 

x + qi 

These constants may be fixed by demanding that the leading asymptotic behavior of the solution 
be unchanged by the backreaction. In particular, for the electric potential we demand the vanishing 
of any 1/x correction to A\ (which would modify the R charge). We also insist that the potential 
continues to vanish at the horizon, A\{xh) = 0, so that remains normalizable at the horizon. 
Note, however, that these two conditions do not fix ki and k2 directly, since /i2 also enters into the 
correction to the potential, as indicated in (j4.3p . Assuming (as we demonstrate below) that /12 falls 
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off faster than 1/x at infinity, the above requirements then lead to 

dx' 



02 (x) 



coth (3i 

12/ ll\ 

V?i (x ) 



— — I dx" x" {x" + qi) X 

{x' + qiy J^, 

2 / // \ 
) 



{xh + qi){x" + Xh + l/g^ + qi)'^ 
From the asymptotic behavior of <pi as given in (j3.20p 

Cl + C2 log X 



h2{xh) 



we obtain 



02 



02(00) + 



coth /?! 



+ 



-C2(ci + C2) + 



+ 



00, 



(4.12) 



(4.13) 



Xh 



xh + qi 



-C2 + 



+ 91 



C2 (Cl + |C2) logX + 



Xh 



■Xh^qx 



4 log2 X + • • • 



(4.14) 



The constant a2{(X)) contributes to a shift in the value of the electric potential (or, equivalently, 
the R charge chemical potential) compared to the zeroth-order solution. 

Turning to the other fields {h2, f2-,X2}-, the coupled system of backreaction equations (j4.5|) . (|4.6|) 
and (|4.8p do not appear to admit a straightforward solution. Nevertheless, some general properties 
are evident. Firstly, in the limit of vanishing electric charge, qi — > 0, the metric equations ()4.5p and 
()4.6p become self-contained. In this limit, (j4.5p may be integrated twice to obtain /12 and, in turn, 
knowledge of /i2 allows (|4.6p to be integrated for /2. Secondly, even with qi 7^ 0, the combination of 
the metric equations admit first integrals, thus reducing the order of the coupled system. In fact, 
using the first-order equation, (j4.9p . along with the backreaction equations for {/12, /2, 02}, we may 
obtain the homogeneous equation 



g^Xh{xh + 1/5^ + qi] 



X 



x + qi 



-h2 



g^Xh{x'^h'2)' + ( x^/2 



„2 W 



-qi coth /?! 



{x + qi){x-Xh) I \ , ^ 2 / , 1 / 

02 H ■ a2 +2g^qi{xh + l/g ) 

Xh + qi 



x + qi 

2^ ^ X — Xh 

x + qi 



X2] =0, (4.15) 



which yields yet another constant of integration (which is presumably related to the conserved 
energy of the static gravitational system) . Finally, we note from (j4.8p that the X2 scalar fluctuations 
satisfy an equation of motion which is of similar form to the <fi equation (jS.Sp . This should not be 
of much surprise, as small scalar fluctuations are naturally governed by the Klein-Gordon equation, 
which in the present background takes on the form 



{x^f^y - \m'^xH^/^^ = 0, 



(4.16) 



for a scalar (Note that, in a background with non-zero R charge, the supergravity potential 
(|2.1ip leads to position-dependent 'masses' for both the (pi and Xi scalars.) Since 99^ and Xi both 
originate from the Tij tensor of the sphere reduction [21], they all have Eq = 2, and hence share a 
common value of mass = —4:g^ (at least when the charges are turned off). Regardless of the 



21 



details, we expect that the X2 fluctuations are hkewise described by a Heun equation. Hence similar 
diflttculties to those we encountered in the previous section arise when obtaining explicit solutions 
for the backreaction of ipi on X2- 

Because the system of backreaction equations is linear, the inhomogeneous solution can in 
principle be obtained by variation of parameters, so long as the fundamental matrix (i.e. the 
complete set of linearly independent solutions to the homogeneous system) is known. More precisely, 
we may rewrite the second-order equations in coupled first-order form. Then the solution to the 
system of n first-order equations 

f'{x) = A{x)x + B{x), (4.17) 



may be formally expressed as 

fix) = $(x)/o + $(x) f ^-^{x')B{x') dx', 

Jxn 



(4.18) 



where $(2;) is the fundamental matrix satisfying $'(x) = A{x)^{x) and normalized according to 
^(xq) = I. This formal solution is not particularly useful, since in the present case #(x) would be 
a 7 X 7 matrix, corresponding to the freedom of four second-order equations ()4.5p through (14. Sh 
minus the first-order constraint (|4.9p . Nevertheless, we note that five out of the seven linearly 
independent solutions to the homogeneous system are easy to obtain: 



1: 

2: 

3: 
4: 
5: 



ho 



f2 



0.2 



X2 



2qi coth Pi 



1 



mx 



X + qi 



X 

(h 

x'^ 

3 
2 



m 

X 

m 
'2^ 



(4.19) 



3m -|- 4qi 



2qix 



1 

X 



-g^{x + qi) + 



7m + Aqi 



2x 



where m = Xh{l + g^{xh + ?)) as obtained from 



The remaining two solutions appear to 



involve non-trivial Heun functions pertaining to the X2 system. 

Our main interest in studying the backreaction is to obtain the asymptotic form of the correc- 
tions so we may discern what effects turning on ipi may have on the conserved quantities such as 
mass and R charge. For the latter, we have already seen from ()4.14p that we may hold the charge 
gi fixed, even as we turn on (pi. For the mass, we turn to an asymptotic expansion of the coupled 
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system (j4.5p . (j4.6p and (j4.8p as x — > oo. From the asymptotic behavior (j4.13p . we obtain 



^9 ~ — 77 



ClC2 - f C2 + g (Xl + ^Xll) 



/2 ~ ^ 



+ (-C2 (ci + 5C2) + giXii) log X - \cl log^ X 



+ 



(xi + \xii) + ic2 (2ci + f C2) 



+ qi 



(ic? + |C1C2 + §cl) 



2 2 2 

3^X11 + 3C2 



2 ^fe(l 
9 



a^h + qi 



C2(3ci + 4c2) logx 



C2 log^ X 



+ 



X2 



(Xi + Xiiloga;) + 



1 



C2 ci + 



Xii 



§£2) + ^ + C2logx 



+ ... 



(4.20) 



Note that here we have used the freedom expressed in (j4.19p to set the leading order terms in /12 
and /2 to zero. The content of the residual two homogeneous solutions are incorporated through 
the constants xi xii; which are related to the shift in the asymptotic profile of the scalar 
Xl at infinity. Of course, the contribution for the full set of homogeneous solutions may have 
to be added back in to satisfy the desired boundary conditions. For the gauge potential, we 
have argued that it is natural to demand A\{xh) = 0. Likewise, here it would be appropriate to 
set f2{xh) + g'^Xhh2{xh) = so that f{x) as given in ()4.2p continues to vanish at x^ when the 
backreaction is included. 

Although we have explored the backreaction of (pi on the other fields of the system, we can 
also see from (j4.20p that the Xi scalar may be deformed as well, through xi and xii- Presumably, 
this would allow a wider class of thermal black hole solutions with scalar hair. However, assuming 
our goal is to thermalize the BPS bubble, one may presumably set xi = Xii = directly without 
any major concern. At the same time, we note an interesting feature of the non-extremal bubble 
solution. While the i?-charged black hole p.4p and the BPS bubble ()2.12p both have the Xi scalar 
satisfying Xi = (for the one-charge case), this condition can no longer be maintained when 

the bubble is thermalized, as evidenced by the non- vanishing backreaction on X2- 

We also note that, according to (j4.20p . it appears that the backreaction can be adjusted so 
that it has no effect on the mass of the black hole. This is because, at least heuristically, the 
mass can be read off from the = 1/x terms in the metric functions Hi and / (since we are 
working in five dimensions), and both such terms are absent in (|4.20p . It may turn out, however, 
that boundary conditions at the horizon would feed in some of the homogeneous solution of the 
third type in (|4.19p . thus leading to a possible shift in the mass. Of course, such issues cannot be 
properly examined in the linearized analysis, and would have to await a full (possibly numerical) 
solution to be addressed. Nevertheless, the heuristic concept of mass can be made rigorous, and 
this is what we turn to next. 
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5 Mass of the non-BPS bubbles 



The definition, and calculation, of the mass of an asymptotically AdS spacetime is somewhat more 
subtle than it is in an asymptotically flat spacetime. One can no longer use the ADM definition, 
which assumes an asymptotically Minkowski region [22]. A generalisation of the ADM procedure, in 
which one decomposes the metric as the sum of an AdS background plus deviations, was introduced 
in [23] by Abbott and Deser. Effectively, one is making an infinite background AdS subtraction 
from a divergent boundary integral. The presence of scalar fields in the solution, such as one has 
in the supergravity black holes and bubbles, can complicate the application of this AD procedure 
considerably, because of the inherent ambiguities in the separation of the metric into background 
plus deviations. Some discussion of the AD approach, and calculations for higher-dimensional 
rotating black holes, can be found in [24]. 

A procedure for calculating the mass of asymptotically AdS spacetimes that avoids all the 
problems inherent in making a split into background and deviation was introduced by Ashtekar, 
Magnon and Das [25,26]. This is based on a conformal definition that expresses the mass in terms 
of an integral of certain components of the Weyl tensor over the spatial conformal boundary at 
infinity. Since the metric approaches AdS asymptotically, the integrand falls off and the integral is 
inherently well-defined. This AMD method was applied in [27] to the calculation of the masses of 
higher-dimensional rotating AdS black holes in general relativity, and in [28] this was extended to 
the case of rotating black holes in gauged supergravities. 

Alternatively, the boundary counterterm method may be used to calculate the mass of config- 
urations in an AdS background [29-33] . This notion of holographic renormalization is particularly 
natural in the context of AdS/CFT, where the addition of boundary counterterms in AdS has a 
natural counterpart in the addition of renormalization counterterms in the dual field theory. Fur- 
thermore, the boundary counterterm method has the advantage that it regulates divergences not 
just in the mass, but also in the on-shell gravitational action which is dual to the thermodynamic 
potential of the CFT. The boundary counterterm method was used in [34-36] to investigate the 
mass of the (undeformed) i?-charged black holes (|2.4p . 

Here, however, we shall use the AMD method to discuss the masses of BPS bubble metrics and 
their non-extremal deformations. Note that, in this definition of mass, pure AdS has by construction 
zero energy. This is in contrast with holographic renormalization, which naturally assigns non-zero 
energy to the vacuum (which is viewed as dual to the Casimir energy of the CFT on x R) . We 
begin by briefiy summarising the AMD procedure, drawing on material presented in [27,28]. 

Consider an asymptotically AdS bulk spacetime {X,g} of dimension D, equipped with a con- 
formal boundary {dX, h}. It admits a conformal compactification {X, g} if X = UdX is the closure 
of X, and the metric g extends smoothly onto X where g = Q'^ g for some function Q with > 
in X and 17 = on dX, with 7^ on dX. One might, for example, take 

n = -, (5.1) 

y 

where in the asymptotic region the metric approaches AdS with R^^ = l^"^ g^u- (For the solutions 
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in gauged supergravity that we shall consider, / = 1/g, where g here denotes the gauge-coupling 
constant.) Since $7 is determined only up to a factor, ^ f Q, where the function / is non-zero on 
dX, the metric g on X and its restriction h = g\dx are defined only up to a non-singular conformal 
factor. The conformal equivalence class {dX,h} is called the conformal boundary of X. If C^upa 
is the Weyl tensor of the conformally rescaled metric g^y = 0,^ g^u, and = d^Q, then in D 
dimensions one defines 

£''y = l'^n^-''^nPn''C^py^. (5.2) 

This is the electric part of the Weyl tensor on the conformal boundary. The conserved charge Q [K] 
associated to the asymptotic Killing vector K is then given by 

where dT,^ is the area element of the (D — 2)-sphere section of the conformal boundary. (The 
derivation of (|5.3p is discussed in [25,26]). Note that the expression (|5.3p is invariant under the 
non-singular conformal transformations of the boundary metric that we discuss above. Thus, one 
may take for any conformal factor that is related to (j5.ip by a non-singular multiplicative factor. 

In order to define the energy, one takes K = d/dt, where t is the time coordinate appearing in 
the asymptotically AdS form 

ds' = -(1 + r') dt^ + -r^^, + dnl_, (5.4) 

of the metric under investigation. The energy (or mass) is then given by 

E = — I E\d%. (5.5) 

Svr (D - 3) Js 



For our present discussion, we need to apply ()5.5p to the class of five-dimensional metrics given 

by 

ds\ = fdt^ + n^l\r^ dr^ + rW^) , (5.6) 

where 7i and / are functions only of r, and dQ"^ is the metric on the unit 3-sphere. It is convenient, 
as usual, to define x = r^. From (|5.5p we then find that the mass is given by 



2^2 



E = lim ^ 



{I- f + xf' -2x^f")n^ + x^{3f'nn' -3fn'^ + 2fnn") , (5.7) 



16/^2 

where primes denote derivatives with respect to x, and we have taken the five dimensional Newton's 
constant G5 = 1. 

Before applying this mass formula to the non-BPS bubbles that we have been investigating in 
this paper, it is instructive to consider some BPS bubble examples. For the spherically-symmetric 
1-charge bubble constructed in [11], the metric functions in ()5.6p are given by 
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where qi is the electric charge parameter. Substituting into (j5.7p . we find that the mass is given by 

E = ^-f. (5.9) 

For the more general case of 3-charge BPS bubbles, which were constructed in [14], we have 

n = HiH2H3, f=l+g^xHiH2H3. (5.10) 

In this case the explicit solution for the three functions Hi is not known, but at large x they take 
the form 

/7i = l + - + •••. (5.11) 

X 

Substituting into (j5.7p . we find that the mass is given by 

E = ^^^S1±^±^, (5.12) 

As a further example, one finds from (j5.7p that the 3-charge non-extremal black holes (j2.4p in 
five-dimensional gauged supergravity have mass given by 

E=^Y.'osh2P, = l(^lm + J2q}j . (5.13) 

In the limit when m is taken to zero, this black hole result reduces to that of the BPS bubble, 
(I5l2]l . 

In general, for a gravitational background parameterized by the asymptotic behavior 

h-i ho 
7Y ^ 1 _L _L J f_ J 



/ ~ l + g^xn + ^ + p^ + -- - , (5.14) 

we find that application of the mass formula (j5.7p gives simply 

E = l {2h, - 3/1) . (5.15) 

This demonstrates that the mass indeed receives contributions only from the 1/x terms in the 
expansion of 7i and /, as alluded to at the end of the previous section. In reality, the expansion in 
()5.14p may also include log terms. However, so long as the logs are confined to the l/x^ and higher 
terms, the mass remains finite and unchanged from (I5.15p . 

For the linearized bubbling AdS black hole solution, the metric functions TC = Hi and / are 
given by the backreaction expansions ()4.2p . As a result, we obtain for the mass 

E = ^ 
4 



3 / 3 

-m + qi+ lim x /i2(x) - -f2{x) 



(5.16) 

Note that we have already assumed that the backreaction functions /12 and /2 vanish at infinity 

/i2(oo) = /2(oo) = 0. (5.17) 
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Based on the backreaction expansion (j4.20p . and taking into account the possible contribution of the 
homogeneous solutions in (j4.19p . we see that the masses of the bubbling AdS black holes remain 
finite, at least for arbitrary linearized deformations. Unfortunately, however, the local analysis 
leading to (j4.20p is insufficient in itself to determine how the mass varies as the ipi deformation is 
turned on while keeping, say, temperature and i?-charge fixed. 

6 A numerical approach 

As in much of the rest of this paper, our interest is in non-BPS bubbles. However, we also numer- 
ically explore some features of the BPS bubbles of [14]. 

6.1 Non-BPS bubbles 

There are two types of AdS bubbles. The first type can be referred to as solitonic AdS bubbles, of 
which the BPS bubbles are a subset. The corresponding geometry is completely regular and horizon- 
free. It interpolates between AdS spacetime at asymptotic infinity and Minkowskian spacetime at 
short distance. The second type is the thermal AdS bubble, for which the geometry contains 
a singular point surrounded by a horizon. In the above sections, we have mainly looked at the 
linearized solution of a subset of thermal AdS bubbles, which we have referred to as bubbling AdS 
black holes. We can demonstrate the existence of both types of AdS bubbles through numerical 
analysis. 



Solitonic AdS bubbles 

Solitonic AdS bubbles have no horizon, and may be described by a global radial coordinate 
X G [0, oo). Focusing on the one-charge spherically symmetric system, the Lagrangian (j2.10p gives 
rise to a coupled set of non-linear ordinary differential equations for the metric functions / and 
Hi as well as the matter fields Xi, ipi and A\. In general, the equations of motion, including 
the Einstein equations, gives rise to five second-order equations (one for each function). However, 
there is also a first-order 'energy' or constraint equation arising from the Einstein equations. As a 
result, any general solution may essentially be specified by nine constants. Not all such solutions 
are physically independent, however. Even with the metric given in the form (|4.ip 

ds" = -H{x)-^/^f{x)dt' + H{xfl^ {-^^ + ^^^i) ' (6.1) 

there remains a residual coordinate transformation 

t^X-^t, x^Xx, H^X'^H, aI^XA], (6.2) 

leaving the form of the metric invariant. This reduces the nine constants to eight physical param- 
eters of the solution. 
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Of course, most of these solutions may be singular as x ^ 0"*". For solitonic AdS bubbles, we 
necessarily demand regularity at the origin, and hence may set up initial conditions by obtaining 
the Taylor expansion of the solution near x ^ 0'^ . For small x, up to linear order, this expansion 
is given by 

Hi = i/o-^^o^o(coshV?-l)^ + --- , 

Xi = Xo + (1 - cosh ip^xl + (cosh^ ifl - l)Xl)x + ■■■ , 

1 2 

cosh ipi = cosh - i(cosh2 - l)H^ {A^ cosh ip'^H^ + 2^X0 - cosh 93? Xq)x H , 



1 



Ai = Ao + ^Ao{cosh^ip'(-1)X^HQ' x + --- , 

_ ^ , (cosh^ ifl - l)AoH^Xo + hI (2 + 4X| cosh (^g + (cosh^ ypg - l)Xg) 

- ^+ 6^0^^ x + ....[b.6) 

Thus, by imposing regularity, we see that the solution ends up being parameterized by only four 
constants, namely Hq, Aq, Xq and (p^. One of the parameters is trivial, owing to the residual 
symmetry (|6.2p . Note that fo = f{x = 0) is fixed to be unity by the equations of motion, and by 
regularity of the spatial slice for shrinking . 

To obtain the BPS bubble, we may impose the first-order condition cosh 991 = (xHi)' given in 
()2.12p . This reduces the four constants down to one according to 

cosh(^? = i7o, ^0 = ^, Xo = H-~\ (6.4) 

and precisely gives rise to the BPS bubbles that preserve ^ supersymmetry, and with charge q = 
Hq — 1. On the other hand, general solitonic (but non-BPS) bubbles may be obtained by relaxing 
any or all of the constraints in ()6.4p . While we do not explore this parameter space in detail, we 
can numerically demonstrate that non-BPS solitonic bubbles can also exist. For example, taking 

Ho = 8, Ao = 1/8, Xo = l/4, cosh(^? = 2 (6.5) 



gives rise to a smooth solution with ipi profile given in Fig. [TJ Taking the scaling symmetry (j6.2p 
into account, these non-BPS bubbles are parameterized by three physical parameters, presumably 
the i?-charge, as well as two scalar 'charges' describing s-wave excitations of the Xi and <pi scalars in 
AdS. These solutions can be viewed as continuous non-BPS deformations of the BPS AdS bubbles. 



Thermal AdS bubbles 

We now turn to AdS bubbles with horizons. Since the horizon appears as a coordinate singu- 
larity, we restrict the numerical solution to cover the outside region x G [xh,oo). In this case, the 
generic regular solution near the horizon may be specified by six parameters, namely the location 
of the horizon x = x^, the initial values of the functions Xi{xh), Hi{xh), cosh ipi{xh) and the slopes 
/i and oi of f{x) and Ai{x), which both much approach zero as x ^ Xh- Note that, taking the 
scaling symmetry (16. 2p into account, we may set Xh to a generic value (say x/i = 1). This indicates 
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Figure 1: Plot of cosh(/?i as a function of x for a non-BPS soliton bubble with parameters given by 
(EH). 



that the solutions may be specified by five physical parameters. However, we find it convenient to 
allow Xh to remain free in the numerical work. 

If we let cosh(/9i = 1, then we turn off the bubble parameter and the solution reduces to that of 
a charged AdS black hole. It is worth remarking that the previously-known charged AdS black hole 
solution, in which the scalar Xi does not have an independent charge parameter, is not the unique 
spherically symmetric black hole. Numerical analysis shows that the generic black hole is, in fact, 
characterized by three parameters: the mass, the i?-charge and the scalar charge of Xi. The thermal 
AdS bubble solution, which also turns on is then characterized by four parameters. (It remains 
unclear what the elusive fifth parameter is. Possibly it could be removed by yet another residual 
gauge transformation.) Using a numerical approach, we can demonstrate the existence of these 
solutions. In particular, we present an AdS bubble solution that can be viewed as a deformation 
of the previously-known AdS black hole. In this case, the initial conditions are specified by the 
following: 

rj . . 1 Y jj-f f l+q + 2x + h 2 Qjxh + 1) „s 

For cosh(/9i = 1, this leads precisely to the previously-known charged AdS black hole. We can turn 
on ipi and numerical analysis indicates that the solution exists provided that 1 < cosh 991 < Hq. In 
Fig. [2] we present a plot of cosh(/7i for q = W, x^ = 1 and cosh = 5. 

6.2 BPS bubbles 

In [14], defining equations for multi-charge BPS bubbles were obtained in D = 4,5,6 and 7 AdS 
gauged supergravities. Those equations in general do not admit analytical solutions. Here we shall 
report that numerical analysis indicates that smooth bubble solutions exist in all of these cases. 
In all of these solutions, the coordinate x (which is denoted as R in [14]) runs from 0, where the 
metric is Minkowskian, to asymptotic AdS spacetime at a; = 00. 
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Figure 2: Plot of cosh(/9i as a function of x for a thermal bubble. We have taken q = 10, Xh = 1 
and cosh ifi = 5. 

In D = 5, the functions Hi describing the general three-charge BPS bubble solutions satisfy the 
nonlinear equation ()2.13p . which may be rewritten as 

where i ^ j ^ k, and where we have defined = xHi. The one-charge case is given by .^2 = = ^) 
while the two-charge case is given by .^3 = x. 

Regularity of the BPS bubble demands that the vanish linearly (or equivalently that the 
Hi approach constants) as a; ^ 0. Taking this into account, we may develop a Taylor expansion 
around x = 0: 

^i = xHi = GiX — ^ajak{af — 1) + ^ai(af — l){8ajak — Oj — a|) + 0{x^) . (6.8) 

To ensure coshc^j > 1, it is necessary that Oj > 1. We verify numerically that smooth solutions 
exist for > 1, and that the functions Hi behave as the following at asymptotic infinity: 

Hi = l + ^-^^ + ... . (6.9) 

Here qi is the charge parameter for the gauge fields Ai , while c\ is the coefficient of the normalizable 
scalar mode according to 

cl + d^logx 

ifi = — \ (6.10) 

X 

at infinity. In particular, the coefficient of the non-normalizable mode always vanishes for these 
BPS bubbles. 

The single-charge case admits an analytic solution of the form given in (|2.15p . This gives the 
exact relation ci = -\/ QiiQi + 2) highlighted in (|2.18p . For the generic multi-charge cases, we have 
been unable to attain an analytical solution to (16. 7p . Nevertheless, we may highlight some of the 
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general features of any such solution. To do so, it is convenient to perform the transformation 
= X + r]i. The resulting differential equation is then 

-1 



x + r]i + 



^x + r]j J \x + 7]k 

The main purpose for introducing this transformation is to arrive at the asymptotic forms 



(6.11) 



r]i ~ (a, — l)x + 0{x^) as x — > 0, ^/i ~ 9« + C ( — ) as x ^ oo. (6-12) 



Note also that it is only the last two terms in (jO.lip which couple the equations for the three fields 
together. 

The AdS vacuum solution is obtained by taking the trivial solution rji = 0. Hence, for small 
deformations of AdS (corresponding to small changes qi <C 1), we may expect r]i ~ 0. More 
precisely, by assuming rn <C x, the above equation may be approximated by 

^ _5M±a. (0.13) 

X + 1 ^ ' 

This approximation also has the feature that the three equations completely decouple in this limit. 
Furthermore, this now admits an exact solution consistent with (j6.12p : 

l+x/(l-Ai) 

m = Ai log 6.14 

1 +x/(l + \i) 

Of course, this is only consistent with our assumption r]i <^ x for Aj <^ 1. In this limit, we obtain 
the approximate solution 

7?,«2A2-^ (A,«l). (6.15) 
1 + X 

Asymptotically, we read off 

q, = 2Xl cl = 2A (Ai«l), (6.16) 

in which case we have demonstrated that 

c\ ~ y% for qi 0. (6.17) 

Physically, we see that since small deformations decouple from each other, we are allowed to turn 
on any independent combination of the three commuting i?-charges {qi, q2 and q^) as we wish, 
while maintaining the form of the 1/8 BPS bubble. 

Larger deformations may be treated numerically. However, it is worth noting the general feature 
of (j6.1ip : since rj starts with a positive slope (forced by demanding coshipi > 1), it will remain 
positive with r]'- < 0. This means r] is monotonic increasing with decreasing slope, and will asymp- 
totically approach its value at infinity, r/j(x — > oo) = qi. The initial slope for r]i, which is — l,then 
determines how large the final charge qi becomes; larger gives larger qi. 
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With this in mind, we may approximate (j6.1ip by noting that the expression + rji) is 
bounded to he between and 1. This expression starts at its minimum value l/a-j when x = and 
increases towards 1 as x — > oo. We now consider the denominator in (|6.11|) : 



den = X + r?, + — — . (6.18) 

\x + r]jj \x + r]kj 

For small x, this is dominated by the initial value of the last term, l/{ajak), while for large x it is 
dominated by the first term. This suggests that we may make the approximation 

den « X + r/j H (6.19) 

in which case the second-order equation 

ri'I = !!Mpl— (6.20) 

x + r]i + [ajak) ^ 

admits a simple solution: 

111 = \Jx^ + 2(^1 + {ajak)~^)x + (ajOfc)-^ - x - {ajakY^ . (6.21) 

This is, in fact, exact for the one-charge case (where we take, e.g. 02 = 03 = 1 and find r/i = 
+ 2{qi + l)x + 1 - X - 1). Here we have chosen the constants of integration in accordance 
with the asymptotic form (j6.12p . The -R-charge is simply gj, and the scalar vevs are 



^i- \|qi{q^ + 2{ajak)-^). (6.22) 
For these expressions to be self-consistent, we examine the behavior of r/j as x — > 0: 

r]i QiUjakX + 0{x^). (6.23) 
Comparing this with (|6.12p . we see that self-consistency demands 

ai = qiajak + l («/j/A;), (6.24) 

which may also be expressed as 

qi = {ai-l)aj^a^\ (6.25) 

where the initial slopes must satisfy a, > 1. This condition gives a range for the allowed values of 
the charges. For example, with three equal charges, the above reduces to 

q = a-^{l - a'^), (6.26) 

which has a maximum value gmax = 1/4, which occurs when a = 1/2. 

Of course, this restriction (|6.25p is entirely contingent on the validity of the approximation 
()6.19p . Numerically, we find that the above captures the qualitative behavior of the solutions, 
but that the actual restrictions on the charges is different. In particular, numerically we find two 
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possibilities for the behavior of qi{ai). For one-charge and two equal charge bubbles, qi can be 
made arbitrarily large. However, for three equal charge bubbles, the actual limiting value of the 
charge is given by 

q < gmax, gmax ~ 0.529. (6.27) 

This cutoff on the maximum amount of ii-charge that can be supported by a three equal charge 
BPS bubble may have implications on the nature of the corresponding 1/8 BPS configuration in 
the dual super Yang- Mills theory. 

For the multi-charge cases where one of the qi can become large, the scalar vev expression (j6.22p 
applies, and we find 

forg^O; 
1 for g ^ oo. 



Qi 



(6.28) 



Even in regions where cannot become arbitrarily large, the dominant behavior is for c\ to 
approach qi from above as qi increases. Similar results can also be obtained for D = 4, 6 and 7 
dimensions, given by 



D 


= 4 : 


Hi 


= Qi- 




D 


= 6 : 


Hi 


= Oj - 


10 v"* 


D 


= 7 : 


Hi 


= Qi- 





(a- - l)aj x2 + iai(a- - l)(4a- - 1) a; + 0{x2) , 



(6.29) 



where the i,j,k indices are not equal when they arise within the same equation. Using this as 
initial data near x = 0, we can find numerically that smooth BPS bubbles exist in all cases for 
O,; > 1. 



7 Conclusions 

We have investigated the non-BPS analog of 'bubbling AdS' geometries in type IIB supergrav- 
ity, corresponding to a special class of non-zero temperature LLM configurations. From the five- 
dimensional point of view, these are solutions of the STU model coupled to three additional scalars 

These solutions can be considered as the bubbling generalizations of non-extremal AdS black 
holes, and have regular horizons. However, unlike the previously-known AdS black holes, the 
bubbling AdS geometries do not have a straightforward generalization away from the BPS limit. 
Thus, we have had to rely on various approximation methods, as well as numerical analysis. In 
particular, we have considered the linearized tpi equations on the background of the non-extremal 
i?-charged AdS black hole. The backreaction onto the other fields occurs at higher order. Even 
at the linearized level, the equations do not have closed-form solutions. Specifically, the linearized 
ifi equations can be mapped to the Heun equation, for which the general two-point connection 
problem remains unsolved. 

Nevertheless, we can find approximate solutions to the Heun equation by matching solutions 
in two overlapping regions. This method is reliable for a certain regime provided that there is a 
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large overlap. In our case, the relevant regime is T S> fii, where T and /x^ are the temperature 
and chemical potentials of the dual thermal Yang-Mills theory, respectively. If, in addition to 
this, we consider a high temperature limit for which T ^ 1, then the ipi equations can be solved 
without matching. Corrections of the order fii/T and 1/T can then be considered via a perturbation 
approach. The complete fi solution is fixed by the boundary conditions at the horizon. Namely, 
we require that (pi is regular at the horizon. For the most part, we focus on the case of one charge 
and a single additional scalar field ipi. 

The behavior of ipi at asymptotic infinity can be related to perturbations away from the UV 
superconformal fixed point of the Yang-Mills theory. In particular, the normalizable mode of ipi 
goes as and corresponds to the dimension two chiral primary operator Tr Zf getting a vev. 
The non-normalizable mode goes as (log r)/r^ and corresponds to a massive deformation of the 
field theory Lagrangian of the form TrZf. We have normalized tpi such that the vev is fixed at 
the value it has for the BPS AdS bubble. Then we discuss how the thermal mass depends on the 
physical parameters of the field theory, namely the temperature, chemical potential and i?-chargc. 
We have considered both the grand canonical ensemble and the canonical ensemble. The fact that 
the Tr Zf term is only present in the Lagrangian for thermal field theories might indicate that there 
is a phase transition at zero temperature. Moreover, we find the thermal mass vanishes at the point 
of the Hawking-Page transition. 

We have gone beyond the linear order in cpi to take into account the backreaction on the other 
fields. In particular, we have used this to obtain the asymptotic form of the corrections in order to 
read off the mass of the bubbling AdS black hole. It would be interesting to further consider the 
backreaction, or else use numerical analysis, to investigate the thermodynamics of bubbling AdS 
black holes. In particular, it would be interesting to find out how the local thermodynamic stability 
constraints, as well as the Hawking-Page transition, may be altered due to (pi. 

It would be useful to understand more concretely whether these bubbling AdS black hole solu- 
tions have any relation to the hyperstar solutions considered in [12]. There, it was suggested that 
the hyperstar background did not have a horizon because the coarse-graining was taken over only 
the half-BPS sector of the full Hilbert space of type IIB theory. Perhaps these bubbling AdS black 
holes are the result of including the non-BPS sector in the coarse-graining. 

Finally, we have shown via numerical analysis that there are actually two types of non-BPS AdS 
bubbles. Thermal AdS bubbles, of which the bubbling AdS black holes are a subset, have an event 
horizon surrounding a singularity. On the other hand, solitonic AdS bubbles are completely regular 
and horizon-free. The latter type of non-BPS bubbles can be obtained from the BPS bubbles by 
smooth deformations, and therefore correspond to non-supersymmetric deformations of the dual 
field theory. It would be interesting to investigate this further. 
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